Practality of Hofstadter Butterfly in Specific Heat Oscillation 
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We calculate thermodynamics! properties of the Hofstadter model using a recently developed 
quantum transfer matrix method. We find intrinsic oscillation features in specific heat that manifest 
the fractal structure of the Hofstadter butterfly We also propose experimental approaches which 
use specific heat as an access to detect the Hofstadter butterfly 
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The interplay between crystalline potential and mag- 
netic field on a two-dimensional electronic gas remains 
a nontrivial problem for decades This issue provides 
a stage on which purely mathematical concepts, i.e., the 
irrationality of a real number, interrupts our intuition of 
physical reality. Hofstadter [2J studied the energy spec- 
trum of the tight-binding limit of this problem, namely, 
the Hofstadter model. He proposed a fractal topology 
for the spectrum(Hofstadter butterfly), which reconciled 
the paradox raised by the irrationality. The experimental 
verification of the Hofstadter butterfly is challenging but 
some hints of the fractal structure have been observed in 
microwave measurements 0, 0|, Hall conductivity [B| and 
magnetic transport measurements Q in analogous sys- 
tems. 

In this paper, we adopt a recently developed quan- 
tum transfer matrix method 7] to study thermodynamic 
properties of the Hofstadter model. We focus on the be- 
havior of internal energy and specific heat as functions 
of magnetic field. As far as we know, this is the first 
report that by theoretical method, the fractal structure 
in the Hofstadter butterfly can be studied by computing 
the specific heat in a magnetic field of a generic value. 
We also briefly discuss the feasibility of experimental ob- 
servations of these features. 

In a previous publication we had used the quantum 
transfer matrix method to study the magnetic proper- 
ties of Hofstadter model. The advantage of this method 
lies in that it directly computes the partition function 
of the model for arbitrary <f>, where <f> is the magnetic 
flux through a unit cell, then the thermodynamic prop- 
erties can be studied steadily. Conventional theoret- 
ical methods, such as Bethe ansatzj^, [l(| and exact 



diagonalizationjlll |12|, are mostly applied to <f> = p/q 



cases, where p and q are mutually prime numbers, and 
q is relatively small. Although detailed information of 
energy spectrum and wavefunction can be obtained with 
these methods, only limited cases of (j) can be studied 
and most discussion was focused on ground state prop- 
erties. Besides, at ground states, due to the fractality 
of the Hofstadter butterfly, the smoothness of physical 
quantities as functions of magnetic field, such as total 



energy, static magnetic susceptibility are significantly di- 
minished. However, within the quantum transfer matrix 
formulation, the effect of finite temperature is embod- 
ied in the partition function at the beginning, and the 
singularities due to the fine fractality will be smeared 
out and the smoothness of physical quantities can be re- 
covered, which makes the comparison to experimental 
results more straightforward. 

Hofstadter model describes the dynamics of two- 
dimensional tight binding electrons in a uniform mag- 
netic field |2|. By applying Landau gauge, i.e., A = 
H(0,x,0), the Hamiltonian is explicitly translationally 
invariant along the y-direction. Fourier transformation 
along the y-axis will then decouple the two-dimensional 
model H into a series summation of one-dimensional 
Hamiltonian H k : 



H = J> 



(1) 



H k = [ tc l,x+l c k,x+tcl !X C k , x+ l 

X 

+2t cos (2nx<f> -k) c\ <x c k , x , (2) 

where k = 2im/N y (n = 0, 1, N y — 1) are the quasimo- 
menta and N y is the lattice dimension along y direction. 
x is the lattice coordinate of electrons along the x-axis. 
<f> is the magnetic flux through each plaquette, with mag- 
netic flux quanta hc/e as unit. 

Hk does not generally have translational invariance 
along the x-axis. But for rational (j) — p/q, periodic- 
ity can be recovered by combining every q cells to form a 
superlattice, and then the problem can be solved by diag- 
onalizing a q x q matrix for each quasimomentum of the 
superlattice. Thus the full energy spectrum and thermo- 
dynamic properties can be steadily obtained, yet appar- 
ently, only up to relatively small q. To study cases with 
a generic </>, the quantum transfer matrix method starts 
from the partition function Z = Tr[exp(— /3H)], which 
can be viewed as an trace of the evolution operator along 
the imaginary time. Since the trace naturally imposes a 
periodical boundary condition, a Fourier transformation 
can be well defined along the imaginary time(the inverse 
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temperature), which is the key point leading us to the 
transfer matrix representation and to significantly sim- 
plify the calculation in Ref.0, H|. Given k, the partition 
function of is defined by 



(3) 



where j3 = The partition function of the whole 

system is simply a product of all ZkS. By making use 
of the translational invariance along the imaginary time, 
Zk can be expressed as a product of N x 2x2 matrices. 
After multiplying different k components, we can obtain 
the partition function of the system, from which one can 
calculate the free energy by F = — (1//3) InZ, and other 
thermodynamic quantities such as magnetic susceptibil- 
ity and specific heat. 

In Ref . |8| , the authors have discussed the effect of lat- 
tice size on the numerical results. Accordingly, we choose 
here N x = 50000 and N y = 100 to ensure the numerical 
accuracy as well as computational efficiency for the tem- 
perature range in this paper. For simplicity, we only con- 
sider the half-filling case, which corresponds to a particle- 
hole symmetry and automatically sets chemical potential 
/i to 0. 

First, we calculate the average internal energy as a 
function of at T = 0.01. As shown in Fig. [IJ at the 
local minima of the internal energy, the electron count 
v{= 0.5 for half-filling) and <6 satisfy the relation in (QJ, 
which was given Ref. ll|. These minima are cusp- like. 



v = M + N<t>, 



M, N e Z. 



(4) 
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FIG. 1: (Color online) The internal energy as a function of <6 
for the Hofstadter model at half filling. T = 0.01. Some values 
of 0(the red number above the top axis) and the integers 
(M, N) corresponding to local minima are marked. 



The global minimum in Fig. Q] is consistent with the con- 
clusion that there is an global minimum of the average 
energy 12j, |l3[ when c/> = v = 1/2, that is, each electron 
carries one flux quanta. We have marked the values of 
c/>(the red number on the top axis) located at distinguish- 
able minima and the corresponding integers M and N in 
Fig. [TJ At zero temperature, the average energy will not 
be smooth almost everywhere because there are infinite 
number of rational cf>s that satisfy (j4|) . But here the tem- 
perature will erase minor singularities and only keep the 
significant ones. 

Then we compute the specific heat from the first order 
derivative of the internal energy with respect to the tem- 
perature. Fig. [2] shows the specific heat C as a function 
of temperature T for some special <fis. The chosen three 
c/>s belong to the pure cases in Hofstadter's proposalQ, 
i.e.,0 = 1/N, or 1 — 1/N when N > 2. Under magnetic 
field of these values, the single Bloch band in zero mag- 
netic field is split into N subbands. If N is odd, the 
central subband has a van-Hove singularity at the center 
point of the energy spectrum(i? = 0). If iV is even, the 
density of states(DOS) goes to zero at E = 0[11|. When 



the temperature is so high that the thermal fluctuations 
are comparable to the energy difference between the low- 
est and the highest subband, the subbands will not be 
able to manifest their internal fine structures from spe- 
cific heat. This can be observed from the high tempera- 
ture tail in Fig. [2] 

The difference in the specific heat for various values 
of 4> will emerge with the decreasing temperature. First, 
at low temperature, the behavior of specific heat can tell 
the singularity of DOS at the energy spectrum center 
point(the Fermi surface(FS) in our half-filling case). In 
the regime near zero temperature, the 0=1/2 and (j> = 




FIG. 2: (Color online) Temperature dependence of the specific 
heat C at half filling for 0=1/2, 1/3, 1/4. 
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1/4 curves are decreasing faster than that of = 1/3. A 
closer observation indicates that = 1/2 and = 1/4 
decrease exponentially-like, while = 1/3 is lincar-likc. 
This is because of the different behavior of DOS at the 
spectrum center point [ll|. For 0=1/2 and = 1/4, the 
original single band in zero field splits up into 2 and 4 
bands. But the centermost two bands are not completely 
separated by a gap, rather they "kiss" at the center point, 
where DOS of both bands goes to zero. Therefore, we 
will expect a gap-like behavior at low temperature, which 
shows up as an exponential-like decrease in specific heat. 
For = 1/3, the Bloch band splits into 3 bands, and DOS 
of the center band is singular at the center point. Thus 
we expect a much slower decrease near zero temperature, 
which is also observed in Fig. [2] 

In the intermediate temperature regime, the specific 
heat tells the information about gaps and redistribution 
of DOS along the energy spectrum. In Fig. [2j the curves 
of = 1/3,1/4 show some similar minor hump struc- 
tures, which is different from the case of = 1/2. For 

= 1/2, two subbands touch at FS, where DOS is zero, 
and there is no finite gap in the energy spectrum. Thus 
there is only one major hump in specific heat curve. For 
both = 1/3 and = 1/4, there is a finite gap lying 
above FS[llJ, which separates the subband on(0 = 1/3) 
or near(0 = 1/4) FS from the higher band, and gives the 
extra minor hump in the specific heat curve. 

Fig. [3r(a) shows the specific heat coefficient (C/T) as 
a function of magnetic field at different temperatures, 
T = 0.1 and T = 0.01. The horizontal axis is chosen as 

1 /0, so that the conventional de Haas- van Alphen(dHvA) 
-like oscillations are shown distinctly in the figure. The 
period A (1/0) of this oscillation is about 2 in both cases, 
which is consistent with that obtained from textbook 
formula [lij]. A(l/0) = A-k 2 /Sf, where Sf is the Fermi 
volume. At half filling, S F = 2tt 2 , thus A(l/0) = 2 in 
Fig.H(a). 

However, a more important observation is that subtle 
oscillations emerge within the dHvA-type period with de- 
creasing temperature and strong field. The specific heat 
at T — 0.1 displays a clean periodic oscillation on the 
weaker field side(large 1/0), while this periodicity is dis- 
turbed in the stronger field regime(small 1/0). This be- 
comes more explicit with lower temperature T = 0.01. In 
the first three periods, very sharp peaks and dips show 
up, and they make a peculiar type of oscillations within 
the period. Even for weaker field regime, some sharp 
structures are still observable. 

For the purpose of comparison, Fig. [3]-(b) shows the 
magnetization oscillation with respect to the magnetic 
field. Similar to the specific heat, the main envelope of 
oscillation is the conventional dHvA oscillation. Besides, 
subtle structures emerge within the dHvA period [8]. By 
comparing the results for magnetization and specific 
heat, we find that the specific heat oscillations are more 
distinct and drastic. 
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FIG. 3: (Color online) (a) The specific heat coefficient C/T as 
a function of 1/0. Two different temperatures T — 0.01,0.1 
are compared; (b) Magnetizationat at T = 0.01 as a function 
of 1/6. 
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FIG. 4: (Color online) Specific heat coefficient (black line) and 
magnetic susceptibility (red line) for the Hofstadter model at 
half-filling. The numerical values of specific heat are 6000 
times larger than the original ones for comparison. The values 
of 0s marked corresponds to some local maxima and minima 
in both specific heat and magnetic susceptibility. 
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To explore the information about the fractal structure 
of Hofstadter butterfly from specific heat, we zoom in the 
first period in Fig J3] and then have Fig[J]for low temper- 
ature specific heat C and magnetic susceptibility x- The 
numerical values of the specific heat are enlarged to 6000 
times of the original values for a clear comparison. Here 
(ft is chosen to be the horizontal axis. The consistency 
between C and x is obvious if we compare the positions 
of local maxima and minima of both quantities. In Fig. [3] 
some fractional values of 0s are marked where they are 
close to the local maxima and minima. Applying Hofs- 
tadter's proposal of constructing the butterfly [2s], we can 
extract the structure of energy spectrum at thess (fts and 
understand why the extrema of C and x are close to 
them. With Hofstadter's proposal, each fractional (ft can 
be decomposed to a set of more "fundamental " fractions, 
or, "local variable" as in Ref.Q, which then directly dis- 
plays the splitting of subbands in the energy spectrum. 
For example, for (ft = 4/13, the center local variable is 
4/5, which means there is a cluster of 5 subbands cen- 
tered at FS. Consequently, a van- Hove singularity of DOS 
shows up at FS and causes the local maximum in C and 
the strong paramagnetism(local maximum in x)0j while 
for (ft = 3/8, the center local variable is 1/2, thus there 
are two subbands lying above and below FS with a zero 
DOS at FS and consequently in FigJH C shows a small 
value(close to 0) around = 3/8 and \ is strongly dia- 
magnetic around (ft = 3/8. 

Therefore, by decreasing the temperature, fractal 
structures of the Hofstadter butterfly manifest them- 
selves by producing peculiar oscillatory features within 
conventional dHvA period. This emergence along with 
decreasing temperature is due to the fact that temper- 
ature provides the only energy scale that sets up the 
resolution of the spectrum. Temperature erases minor 
bands and gaps that are smaller than the scale of temper- 
ature and restores the smoothness of physical quantities. 
But fractal structures with an energy scale larger than 
the temperature survive, and are able to manifest them- 
selves by displaying smoothened singularities in thermo- 
dynamic quantities. Thus the subtler fractal structures of 
Hofstadter butterfly can be probed by the measurement 
of the specific heat at lower temperatures. 

We propose to adopt the superconducting thin 
films(for example the element Nb) with periodic arrays 
of pinning sites [15| to realize this temperature-dependent 
emergence of fractal structures in specific heat. The ar- 
tificial pinning centers hold great potential. Just below 
the onset temperature of superconducting transition, the 
electrons possess long mean free path. When the interval 
between adjacent sites comes to the order of lOOnm, the 
experimentally accessible steady fields can enter the in- 
teresting regime of (ft. The resulting effective lattice sub- 
jected to perpendicular magnetic field is probably able to 
show the fractal properties of the Hofstadter model. In 



addition, the purity requirement of the sample is relaxed 
when considering the specific heat measurement. 

In summary, adopting the quantum transfer matrix 
method, we compute the internal energy and specific heat 
of the Hofstadter model, and for the first time we study 
the oscillation of the specific heat with varying magnetic 
field as a signature of fractal structure of the Hofstadter 
butterfly. In low field regime, the oscillation period of 
specific heat is consistent with the conventional dHvA 
oscillation. When the temperature is decreased, sharp 
peaks and dips emerge in addition to the dHvA-type 
background. These peculiar oscillatory behaviors are di- 
rect indications of DOS in the fractal energy spectrum. 
We also suggest the possibility of making use of supercon- 
ducting films to detect this fractal structure by measuring 
the specific heat. 
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